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We show that there exists a systematic expansion around four spatial dimensions for Fermi gas 
in the unitarity regime. We perform the calculations to leading and next-to-leading orders in the 
expansion over e = 4 — d, where d is the dimensionality of space. We find the ratio of chemical 
potential and Fermi energy to be fi/ep =2^ +16^ Ine — 0.0246 e ' + ■ ■ ■ and the ratio of the 
gap in the fermion quasiparticle spectrum and the chemical potential to be A//i = 2e'^ — 0.691 + - ■ • • 
The minimum of the fermion dispersion curve is located at |p| = (2meo)^ where eo/m = 2 -|- 0(e). 
Extrapolation to d = 3 gives results consistent with Monte Carlo simulations. 
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Introduction. — Dilute Fermi gas at infinite scattering 
length 0,|3| has attracted considerable attention recently. 
The system can be realized in atomic traps using the 
Feshbach resonance |1I1I1I1[1I1I3- It might be rel- 
evant for the physics of dilute neutron gas [lO|. It has 
been suggested that its understanding may be important 
for the understanding of high-Tc superconductivity [111 . 
From the theoretical perspective this is a unique nonrel- 
ativistic system that has no intrinsic scale parameter. 

Theoretical treatment of the system is difficult, how- 
ever, precisely due to the lack of the any small dimension- 
less quantity. The usual Green's function techniques of 
many-body physics become completely unreliable since 
the expansion parameter is large, na^ ^1. So far, no 
systematic treatment has emerged. Recently consider- 
able prog ress has been made by Monte Carlo simula- 
tions [ij, llJj llJi llH ■ However, there are many reasons 
that make an analytical treatment, if it exists, extremely 
useful. First, there are many problems that still can- 
not be solved by Monte Carlo simulations. Examples 
include polarized Fermi gases, recently realized in exper- 
iments [l6|, ll7|, but whose lattice realization suffers a 
fermion sign problem, and questions related to dynamics 
like the dynamical response function and the kinetic co- 
efficients. Second, in many cases analytical approaches 
give unique insights that are not obvious from numerics. 

In this Letter we propose an approach based on an 
expansion around four spatial dimensions. In this ap- 
proach, one would be doing calculation in 4 — e spatial 
dimensions, where the small number e is used as a pa- 
rameter of the perturbative expansion. Results for the 
physical case of three spatial dimensions are obtained by 
extrapolating the series expansions to e = 1. This ap- 
proach has been extremely fruitful in the theory of the 
second order phase transition [l8j. In our case, we find 
that even at e = 1 the series over e is reasonably well- 
behaved, strongly suggesting that the limit d ^- 4 is not 
only theoretically interesting but also practically useful. 

The special role of four spatial dimensions has been 
recognized by Nussinov and Nussinov 1W\. They noticed 



that at infinite scattering length, the two-body wavefunc- 
tion has a 1/r^ behavior when the separation between 
two fermions r becomes small. The normalization inte- 
gral of the wavefunction has a logarithmic singularity at 
r ^ 0, from which it is concluded that at d = 4 the sys- 
tem must become a noninter acting Bose gas. As far as we 
know, no other attempt to exploit this special property 
of four dimensions has been made prior to our work. 

Feynman rules and the counting of the powers of e. — 
Due to the universality of the unitary Fermi gas, any 
short-range two-body interaction can be used, if it cor- 
responds to the infinite scattering length. In particular, 
we can choose to work with the Lagrangian of local four- 
Fermi interaction. After a Hubbard-Stratonovich trans- 
formation, the Lagrangian density of the unitary Fermi 
gas can be written as (here and below h = 1) 



/: = *t 



idt + ^ — * + M*V3* 

Co 



(1) 



where cq is chosen to correspond to infinite scattering 
length. In dimensional regularization, which we will use, 
Co = 00. From now on we set 1/co = 0. Here ^ is 
a two-component Nambu-Gor'kov field, "f = (V't,V'|)'^j 
f± — ^('^1 ^ i'^2), and 171.2,3 a-re the Pauli matrices. 

The ground state is a superfiuid state where (j) con- 
denses: (0) — 4>o. We choose (j>o to be real. Then we 
expand 



cj) = cj)o + gif, g = H;— ' 2 

m \ ZTT J 



where g ~ 0(e^' ^) was chosen for later convenience, and 
rewrite the Lagrangian density as a sum: C = Cq -{- Ci -\- 
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FIG. 1: Feynman rules. The two vertices on the last column 
come from £2, while the rest from C\. Solid (dotted) lines 
represent the fermion (boson) propagator iG (iD). 



C2, where 

£0 = *^ (idt + ^ + O-+0O + CT_0o ) * 

V 2m 



+ '^*(^^* + ^)'^' 



(3) 



Ci = 5r*V+*(^ + .g*V_\I/(^* + fi'^'^as'^ + 2fiip*ip, (4) 



(5) 



£2 = -V^* ( idt + -7— ] f - 2^tp*if . 



The part Cq is a Lagrangian density of noninteracting 
fermion quasiparticles and a boson with mass 2to, whose 
kinetic terms are introduced by hand in Co and taken 
out in £2- The propagators are generated by Co and the 
vertices by Ci and C2 ■ The fermion propagator is a 2 x 2 
matrix, 



G{po,p) 



pI- El + i5 \ -00 Po - £p 



(6) 



where £p = p^ /2m and Ep = [e. 
propagator is 



t>lY'^- 



D{po..p) 



Pi) 



2 



i5 



The boson 



(7) 



The vertices come from Li and C2 and are depicted in 
Fig. n where 



no=Po-f. 



(8) 



The fermion-boson couphng is proportional to g and is 
small in the limit e ^ 0. 

Let us first consider Feynman diagrams constructed 
from Co and Ci only, without the vertices from £2- We 
make a prior assumption /i/0o ^ e, which will be checked, 
and consider 0o to be 0(1). Each pair of boson- fermion 
vertices brings a factor of e, as each /x insertion. Therefore 
the naive power of e for a given diagram is Ng/2 + N^, 
where Ng is the number of vertices and N^ is the number 
of /i insertions. However, this naive counting does not 
take into account the fact that there might be inverse 
powers of e coming from integrals which diverge at d = 4. 
Using a power counting similar to that in relativistic field 



theories, one can show that inverse powers of e appear 
only in diagrams with no more than three external legs. 
Moreover, from the analytic properties of G{p) and D{p) 
in the ultraviolet region, one can show that there are 
only four diagrams which have 1/e singularity near four 
dimensions. They are one-loop diagrams of the boson 
self-energy [Figs. |2Ia) and|2Ic)], ip tadpole [Fig. [Sfe)], 
and vacuum (the middle of Fig. ^. The diagrams in 
Figs.|2Ia) and|2Ic) combine with the vertices from £2 to 
restore the naive e power counting. 

For example, the diagram in Fig. |2Ia) is 



''^»^'-''Jwk^°"[>'-'2)H 



^-i 



(9) 

The integral has a pole at d = 4, so it is 0(1) instead of 
0(e) according to the naive counting. The residue at the 
pole can be computed as 



n(p) 



Po 



Oie), 



(10) 



which is canceled out exactly by the vertex IIo in £2. 
Therefore the sum of Figs.|2la) andEl^b) is 0(e). 

Similarly, the diagram in Fig. El^c) contains a 1/e sin- 
gularity, and is 0(e) instead of naive O(e^). The leading 
part of this diagram is canceled out by the second vertex 
from £2, and the total is again O(e^). 

Finally, the ip tadpole diagram with one fi insertion 
[Fig.[2i;e)] is 0(ei/2) instead of naive ©(e^/^). The only 
diagram that can cancel this is the tadpole diagram with 
no p, insertion. Fig. I3f). The condition of cancellation 
determines (poip) to leading order in e. This condition 
will be automatically satisfied by the minimization of the 
effective potential. 

Thus, we can now develop a diagrammatic technique 
for our system. For any Green's function, we write down 
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FIG. 2: Restoration of naive e counting for the boson self- 
energy and the cancellation of tadpole diagrams. The fermion 
loop in (c) goes around clockwise and counterclockwise. 
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FIG. 4: One- loop diagrams for the fermion self-energy of order 
0{e). 



FIG. 3: Vacuum diagrams for the effective potential up to 
next-to-leading order in e. The second diagram is 0(1) in- 
stead of naive 0{e) because of the 1/e singularity. 



all Feynman diagrams according to the Feynman rules, 
using the propagators from Cq and the vertices from Ci. 
If there is any subdiagram of the type in Fig. |2Ia) and 
Fig. [2Ic) , we add a diagram with a vertex from C2 ■ The 
result will be 0(6^9/2+^^) [2l|. 

Leading and next-to-leading order results. — We shall 
perform explicit calculations, employing the Feynman 
rules and the e power counting that have just been devel- 
oped, to leading and next-to-leading orders. The depen- 
dence of (f>o on /i is most conveniently computed from the 
minimization of the effective potential Vcff{4>o) |20j. To 
next-to-leading order, the effective potential receives con- 
tribution from three vacuum diagrams drawn in Fig. |31 
fermion loops with and without a /i insertion and a 
fermion loop with the boson exchange. The contribution 
from the one-loop diagrams reads 



^i(0o) = ^ 



[. , 7- 


-3(7 + ln2) " 


f m(f>Q 


i-H 


6 'J 

-2(7-ln2) " 
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Iv 27r 



d/2 



d/2 



(11) 



where 7 « 0.57722 is the Euler-Mascheroni constant. 
The contribution of the two-loop diagram is 



V-2 



2l<P0 



)=9' 



^^%^G,MG22{q)D{p-q) (12) 



(27r 
g'^ f dpdq {Ep - ep){Eq ~ Sq) 



4 J {2TT)^d EpEq{Ep + Eq+ £p-q/2) ' 

This integral is convergent even at d = 4. Its value is 

d/2 



V2{4>o) = -Ce 



2tt 



(13) 



where the constant C is given by a two-dimensional in- 
tegral 



C 



dx / dy 
Jo 



[f{x)-x][f{v)-v] 

fi^)f{y) 



X 9ix,y) - \/g^{x,y) - 



xy 



(14) 



with /(x) = (x2 + l)i/2 and 5(x,7;) = /(x)-f/(2/) + i(.x + 
y). The result of the numerical integration is 



C « 0.14424. 



(15) 



The minimum of the effective potential Vcsift^o) 
ViiM + ^2(^0) is located at 



2^ 



[l + (3C-l-f ln2)e]. 



(16) 



Note that the previously made assumption fi/4>o — 0(e) 
is now checked. Also if one used the mean field approx- 
imation, one would reproduce the leading 2/i/e term in 
Eq. H16|l . but not the 0(e) correction. The value of Vcs 
at 00 in Eq- l|16|l determines the pressure P = ^Kff (0o) 
at chemical potential /i. The density is determined from 
n — dP/d^, and the Fermi energy from the thermody- 
namic of free gas in d dimensions is given by 



27r 
Ef = — 

771 



1 fd 
2M2 



2/d 



l-ln2 



The nontrivial power of e comes from taking 71 ~ e 
the power of 2/d. We find the parameter ^ = n/e-p, 



(17) 
' to 



, e3/2 / elne 



1-(3C--(1-In2))6 



Substituting the numerical value for C, one finds 



C 



1 



3/2 



\:t''^\n€-Qmme^'^ + --- 



16 



(18) 



(19) 



The smallness of the coefficient in front of e^/^ jg ^ re- 
sult of a cancellation between the two- loop correction and 
the subleading terms from the expansion of the one-loop 
diagrams around d = 4. 

Quasiparticle spectrum. — To leading order, the dis- 
persion relation of fermion quasiparticles is w(p) = Ep = 
(ep-|-(/)o)^^^. It has a minimum at p = with a gap equal 
to A = 00 = 2/i/e. The next-to-leading order correction 
comes from three sources: from the correction of 0o in 
Eq. I|16(l . from a p, insertion to the fermion propagator, 
and from the one-loop self-energy diagrams, —iY,{p), de- 
picted in Fig. 0] Using the Feynman rules one can see 
that there are corrections only to the diagonal elements 
of the self energy: 



Sii(p) = ig 



dk 



■G22{k)D{p-k) 
Ek 



{2ttY+^ 
(f_ r _dk 

'2 J {2ttY Ek{Ek+ek-p/2-po) 



£fc 



(20) 



and T,22iPo,p) = ~^iii~Po,p)- To find the correction to 
the dispersion relation around its minimum, we only have 



to evaluate the self-energy atpo = 4>o,ep^fi and expand 
E(po,p) = SO(0o,O) + I]'((^o,0)£p/(/.o. By solving the 
equation det[G'~^(w,p) + ficr^ — E(a;,p)] = in terms of 
to, we see the dispersion relation around its minimum is 
given by 



uj{p) = A 



sof 



200 



(21) 



where A ^ <j>q + (Ej^ + S^2)/2 and eq = M + (S^2 - 
E0i)/2-(E'ii + E^2)/2. 

The result of an explicit calculation is the following. 
The minimum of the dispersion curve is located at a 
nonzero value of momentum, \p\ = (2?7ieo)^ j where 



£0 = 2/x. 



(22) 



Note the difference with the mean field approximation, 
in which Eq — fJ-- The correction to the gap is 

^[Sn(0o,O) + S22(<^o,O)] = -e(81n3-121n2)(/.o. 

(23) 
Combining it with the correction in Eq. (|16|l , we obtain 






(3C- 1 -81n3 + 131n2)e] 



- - 0.691. 

. (^^) 
Extrapolation to e — 1. — Although the formalism is 

based on the smallness of e, we see that even at e = 1 the 
corrections are reasonably small. If we naively use only 
the leading and next-to-leading order results, extrapola- 
tion to e = 1 gives for three spatial dimensions 



C « 0.475, 



£0 ^ 


s 2 


i\ 


M 




M 



1.31. 



(25) 



They are reasonably close to the results of recent Monte 
Carlo simulations, which yield £, ~ 0.42, £o/m ~ 1-9j 
and A//i « 1.2 |15||. They are also consistent with recent 
experimental measurements off, where f = 0.51±0.04 9| 
and f = 0.46 ±0.05 [l2|. Thus there is a strong indication 
that the e expansion is useful in practice. A calculation 
of the e^ corrections to these results would be extremely 
interesting. 

Conclusion. — We have developed a systematic expan- 
sion, treating the dimensionality of space as close to 
four, and obtained very reasonable results. As far as 
we know, this is the only systematic expansion for the 
unitary Fermi gas at zero temperature that exists at this 
moment. We found that the the corrections are not too 
big even when extrapolated to e — 1, which suggests 
that the picture of the unitary Fermi gas as a collection 
of weakly interacting fermionic and bosonic quasiparti- 
cles may be a useful starting point even in three spatial 



dimensions. There is a host of problems that can be ad- 
dressed using this approach: the phase diagram of the 
polarized system, the structure of the superfluid vortex, 
finite-temperature physics, etc. It is interesting to note 
that the critical dimension of a superfluid-normal phase 
transition is also four, making weak-coupling calculations 
reliable at any temperature for small e. 
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son propagator is needed to avoid infrared singularities. 
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